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Propagation  Study  for  a  Tropospheric  Transhorizon 
Radar 

KURT  TOMAN,  SENIOR  MEMBER,  IEEE 

Derived  from  an  early  National  Bureau  of  Standards  (NBS) 
tropospheric  transhorizon  propagation  data  base  [1],  a  par¬ 
ticular  loss  term  defined  as  path  attenuation  was  used  in  the 
radar  equation  to  estimate  the  behavior  of  signal-to-noise  ratio 
with  frequency  (10  to  1000  MHz)  and  distance  (SO  to  1000  km) 
including  median  and  1,  10,  90,  99  percent  variability.  This  loss 
term  depends  on  frequency,  distance,  climate,  refractivity  and 
an  empirically  derived  attenuation  function  which  takes  into 
account  effects  of  antenna  heights.  Equivalent  antenna  temper¬ 
atures  are  assumed  to  be  due  to  galactic  noise.  Receiver  noise  is 
characterized  by  noise  figure.  For  conveniently  chosen  system 
parameters  (1  MW  transmitter  power,  1  Hz  receiver  bandwidth, 
antenna  gains  G,  =  Gr  =  1000),  signal-to-noise  ratio  behavior  is 
illustrated  for  several  target  cross  sections  (e.g.  resonant  dipole, 
conducting  sphere)  as  a  function  of  frequency,  distance,  varia¬ 
bility  for  particular  climates.  Conversion  to  other  system  param¬ 
eter  values  is  straightforward.  The  complete  details  are  available 
from  the  author  (2] . 
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Optimization  of  Radio  Communication  in  Media  with 
Three  Layers 

GERVASIO  P.  S.  CAVALCANTE  AND  ATTILIO  J.  GIAROLA. 

SENIOR  MEMBER,  IEEE 

Abstract— Tht  electromagnetic  radiation  of  an  electric  dipole  in  a 
medium  with  three  layers  is  examined  using  dyadic  (ifeen’s  functions. 
The  far  zone  field  for  problems  of  this  nature  is  primarily  determined 
from  the  lateral  wave.  It  is  shown  that  the  excitation  of  this  wave  may  be 
reinforced  through  a  dipole  inclination  and  an  optimum  position  may  be 
determined.  The  radio  losses  for  typical  forests  were  calculated  for 
vertical  and  horizontal  dipoles  and  for  dipoles  with  an  optimum  in¬ 
clination.  The  theoretical  results  are  in  good  agreement  with  the  available 
experimental  data. 


1.  INTRODUCTION 


z 


The  analysis  of  the  electromagnetic  wave  propagation  in 
layered  media  has  been  considered  by  many  authors  [  1 )  -[7] .  A 
model  commonly  used  to  analyze  the  radio  propagation  in 
forests  consists  of  a  homogeneous  and  isotropic  dielectric  layer 
placed  over  a  conducting  flat  earth  [8] -[  10] .  In  this  analysis 
the  Hertz  potential  theory  is  used,  and  particular  attention  is 
devoted  to  the  behavior  of  the  lateral  wave. 

The  use  of  the  dyadic  Green’s  functions  for  apnlysis  of 
the  electromagnetic  wave  propagation  in  semi-,  j.  media 
was  described  by  Tai  [11]  and  a  generalization  o  diese  func¬ 
tions  for  the  case  of  an  AMayered  medium  has  beep  recently  ob¬ 
tained  [12]. 

In  this  communication  the  dyadic  Green  *  functions  for  a 
three-layered  medium  are  determined.  From  these  functions  the 
electric  fields  of  the  electromagnetic  waves  radiated  from  the 
vertical  and  horizontal  electric  dipoles  were  obtained.  In  addi¬ 
tion,  the  electric  field  of  the  electromagnetic  wave  radiated  from 
an  inclined  electric  dipole  was  obtained,  and  an  analysis  was  de¬ 
veloped  to  evaluate  the  inclination  angle  for  optimum  excitation 
of  the  lateral  wave.  The  results  are  expressed  in  terms  of  radio 
losses  between  transmitter  and  receiver  and  for  typical  forests. 

II.  OUTLINE  OF  THE  GENERAL  THEORY 

The  geometry  that  will  be  considered  here  is  shown  in  Fig.  1 . 
An  electric  dipole  with  an  inclination  a  with  respect  to  the 
boundary  plane  between  layers  and  an  observation  point  are  both 
located  in  the  intermediate  layer  (medium  2).  The  height  of  this 
layer  measured  along  the  z-direction  is  H,  while  the  height  of  the 
other  two  adjacent  layers  is  infinite. 

The  source  that  feeds  the  electric  dipole  has  a  harmonic  time 
dependence  given  by  exp  (-/w r).  The  dielectric  permittivities  of 
the  three  media  are  expressed  in  their  complex  form,  e„  = 
CrnC 0  =  +  Mn  =  e«(  1  +  jo„/coe'n),  where  ern  is  the  complex 
dielectric  constant  and  o„  is  the  conductivity  of  medium  rt(n  = 
1, 2,or3),e0  is  the  free-space  permittivity.  The  real  and  imaginary 
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Fig.  1.  Geometry  for  model  with  three-layered  medium. 

parts  of  the  dielectric  permittivity  of  medium  n  are  given  by  e', 
and  respectively,  and  the  wavenumber  is  given  by  k„  = 
(jjs/no^n-  The  magnetic  permeability  of  any  one  of  the  three 
media  is  assumed  to  be  equal  to  the  free-space  permeability  n0. 

The  wave  equation  for  the  electric  field  Es  in  layer  s  due  to  a 
source  in  layer  f(s,f=  1 , 2,  or  3)  is  given  by  (see  Fig.  1): 

VxVxE,  -kjE,=  juHoJfSfs,  ( 1 ) 

where  is  the  Kroneker  delta  (6{  equal  to  one  for  f  =  s  and 
equal  to  zero  for  /  ¥=  s)  and  Jf  is  the  current  density  from  the 
source. 

The  electric  field  at  a  position  R  from  a  source  at  R  may  be 
obtained  from  a  dyadic  Green’s  function  of  the  electric  type, 
G^\RIR'),  that  satisfies  the  following  wave  equation; 

VxVx8isn(R/R')  -  k 2Giit\R/R')  =  /6(R  -  R'yb{.  (2) 

where  16  (R  -  R')  is  the  unit  delta  dyadic  function  and  the  su¬ 
perscripts  ( sf)  refer  to  the  layer  where  the  observation  point  and 
the  source  position  allocated,  respectively. 

Knowing  G^\R/R')  and  J/(R'),  the  electric  field  may  be  ob¬ 
tained  from  (see  Tai  [11],  for  example); 

E^R)=itJH0  Iff  8isf)(R/R')-jf(R')dV  (3) 

*7 

where  the  integral  is  over  the  whole  volume  Vf  that  contains  the 
source  in  layer/and  G^\R/R')  satisfies  the  following  boundary 
conditions  (for  convenience,  the  functional  dependence  (/?//?') 
will  be  omitted  in  all  dyadic  Green’s  functions); 

zx8<sf)  =  zx8(/*l'n  (4) 

zxMxGW  =  zxVx8</+ 1 S*  (5) 

at  the  interfaces  z  =  0  and  z  =  H.  in  addition,  G^n  satisfies  the 
Sommerfeld’s  radiation  condition. 

The  method  of  scattering  superposition  (see  Tai  [11],  for  ex¬ 
ample)  is  used  to  obtain  G^\  Thus 

5^  =  +  G(/n.  (6) 

The  dyadic  Green’s  functions  may  be  expressed  in  terms  of 
the  cylindrical  vector  wave  functions  MCnX(h)  and  NtnK(h)  (see 
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Tai  [11],  for  example),  where  the  subscripts  e  and  o  indicate 
even  and  odd  functions,  respectively,  n,  X,  and  h  are  the  eigen¬ 
values  associated  with  the  0,  r,  and  z  coordinates,  respectively. 
The  eigenvalues  h  and  X  are  related  by  It  =  \Jk 2  -  XJ . 

For  simplicity,  the  subscripts  of  the  vector  wave  functions 
Menk(h)  and  NtnK(h)  will  be  omitted.  The  prime  on  these  func¬ 
tions  M(h),  N  (It)  is  used  to  indicate  that  they  are  expressed  in 
terms  of  the  coordinates  (r ,  0',  z'). 

The  free-space  dyadic  Green’s  function  is  given  by  [13] 


G0 


zzb(R  “*')+/  J- 
ki  4jr 


’  S  (2-«o) 

n=0 


M(h2)Af(-h2)  +N(h2)N,(-h2),  z>z 
M(-h2yft(h2)  +  N(-h2yf(h2),  a  <  z  \ 


From  (3),  „  I),  and  (11)  one  obtains,  for  the  region  H>z>z': 

E2(R)  =  -^  j0  £[px|[AW(*2X1+«) 

+  cM„ ,  x(-*2)]  e,h2Z'  +  W Mo  i  x(h2) 

+  aMolK(-h2)]e-fh2 ■*' 

+  [**,a(*2XI  +  b)  +  dNelK(-h2)] 


(-/  £  V"2*'  +  [»'*.■  xfc)  +  bNelK(-h2)  | 

(/  J;)*"7*2*  j  oM 

[(1  +  b)e/hix'  +  b'e-lhit] 


(7) 

where  50  is  the  Kroneker  delta  (60  =  1  for  n  =  0  and  60  =  0 
for  n  #  0).  Since  the  interest  of  this  communication  is  for  fields 
away  from  the  source  region  the  first  term  on  the  right  side  of 
(7)  is  irrelevant  and  will  be  omitted  from  the  expressions  of  the 
dyadic  Green’s  functions. 

The  scattering  dyadic  Green’s  functions  are  conven¬ 

iently  chosen  so  that  the  boundary  and  radiation  conditions  for 
^  achievable.  The  result  is 

^r12)  =  7-  I  2  (2-60){M(fi,)[eAf'(-A2)  + 

47T  #q  /I2A  n  =  o 

+  e'Af(h2)]  +  N(h  j  )(fff(-h2)  +  /'Atyi2)] }  (8) 

d\  *** 

6<22>  =  -/-  /  —  2  (2-«o){^2)[aW'(-*2) 

4tt  /q  /IjA  n=o 

+  a'F(/i2)]  +  N(h2)[brf(-h2)  +  b'N'ih^] 

+  M(-h2)[cM\-h2)  +  aM?(h2)] 

+  N(-h2)(dN'(-h2)  +  bN(h2)\ }  (9) 


^”>  =  4"  /  fT  2  (2  -«o) 

47T  /q  /I2A  n-Q 

•  {A/(-/<3)[fr'^2)  +  «M'(-*2)] 
+  N(-h3[h,N‘(-h2)  +  hN'(h2)\} 


+  A.o  x(-A2)[*/'’2 


*'  +  f»e-/',2*']} 


(12) 


where 


exp  (/2A2//)  ,  R? 


exp(/2/i2tf) 


6  = 


,,  R?  exp  (j2h2H)  ,  R ^  exp  (/2/i2//) 

0  C  ft  u  |» 

£>^  dP  Dv 


03) 


with 


h2—h  1  r>H  h2—h3  nV  k\h2—k\hi 

A  1  —  A!  —  A  |  —  I  I 

A2+fc,  h2  +  fr3  k}h2+k\hx 

R2=lfuZfu  D^'V—  1  —  Rlx-yR2'V  exp  (j2h2H). 
k\h2  +  k\h3 

(14) 

Using  the  expressions  for  MtnX(h)  and  NtnX(h)  ( see  Tai  [11], 
for  example)  the  z  component  of  (12)  may  be  obtained  as 
1  +R\el2»2(H-^ 


r/XX2 


(10) 


where  the  coefficients  that  appear  in  (8)-<  1 0)  are  obtained 
from  boundary  conditions  (4)  and  (5). 

The  current  density  of  an  electric  dipole  with  an  inclina¬ 
tion  a  with  respect  to  the  boundary  plane  between  layers,  as 
shown  in  Fig.  1 ,  may  be  expressed  as 

j 2(R')  =  (Pxx  +  pzz)[6(x  -  0)60-  -  0)6(z  -  z’)]  (11) 

where  px  and  pz  are  the  horizontal  and  vertical  dipole  moments, 
respectively. 


#r  _  ■  ^0  f 

r_,4 nk\  J0  Unn  1  -R\R\e'2^H 

’  |  ^Px  cos  07 1  (Xr)  +  jpz  ■—  70(Xr)J  [1  +R%e/2h*H] 

05) 

III.  z -COMPONENT  OF  THE  RADIATED  ELECTRIC  FIELD 

The  integrand  of  (IS)  has  two  terms,  one  proportional  to  the 
vertical  dipole  pt  and  the  other  to  the  horizontal  dipole  px.  Each 
integral  will  be  evaluated  separately.  The  z-component  of  the 
radiated  electric  field  due  to  the  vertical  dipole  is  obtained  from 

/“  X3 

0(X)-  wy>(Xr)rfX  06) 
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where 
Q  = 


“PoPz 

Snk2 


(17) 


m= 


[e~/h2*  +R%elh*z  ] ^lhi*  + 

[1  -R\R\e‘lh2H] 


(18) 


and  the  Bessel  function  JqQj)  was  transformed  into  a  Hankel 
function  H^(\r)  in  order  to  change  the  limits  of  integration 
from  zero  to  infinity  to  the  more  convenient  limits  from  minus 
infinity  to  infinity. 

The  integral  (16)  represents  the  field  in  medium  2  and  results 
from  the  superposition  of  cylindrical  waves  thar  are  successively 
reflected  at  the  interfaces  between  media  1  and  2,  and  media  2 
and  3,  with  reflection  coefficients  R ,  and  R2 ,  respectively. 

Since  the  interest  here  is  the  calculation  of  the  fields  in  the  far 
zone,  the  Hankel  function  which  appears  in  (16)  may  be  replaced 
by  its  asymptotic  expression.  In  addition,  one  over  the  expression 
that  appears  in  the  denominator  of  (18)  may  be  expressed  in  an 
infinite  series  form  (binomial  expansion).  This  however,  does  not 
remove  the  necessity  of  considering  the  branch  points,  as  will  be 
done  afterwards.  After  exchanging  the  integration  and  the  sum¬ 
mation  operations  and  using  the  transformation  X  =  k2  sin  v,  (16) 
may  be  written  as 


Ez  =  2  <t'‘m  09) 

m-0 

where 

lm=  /  {e-b'  +R!el>W  +  R\el’<2H-") 

•  (R\  R2e2bH)m  •  sin3/J  v  •  elk2riin"  dv,  (20) 

(21) 

b  =  jk2  cos  v.  (22) 


field  due  to  a  horizontal  polarized  source  may  be  obtained.  The 
superposition  of  both  fields  leads  to  the  expression  of  the  z 
component  of  the  electric  field  due  to  an  inclined  electric  dipole. 
The  result  is 


Ei=-Et  2  [/?2>i)e“,/]mll  +/n(l  +  *-«<«-*>)] 

m  =  0 


[l+/?2>!)e 


r  «1  COS0COSO 

■***]•  - sin  a  H - — - - 

[_1 —  «i  Vl  —  n2 


where 


VE[  +  P J 


pz  sin  a 


vr 

(V/m)  (27) 

(28) 


IV.  OPTIMUM  INCLINATION  ANGLE 


As  one  has  concluded  previously,  the  most  important  contri¬ 
bution  to  the  radiated  field  results  from  the  lateral  wave.  An  im¬ 
portant  characteristic  of  this  wave  is  that  its  excitation  is  asso¬ 
ciated  with  the  angle  of  incidence  for  total  reflection,  vc  = 
sin-1  (ri|).  Thus  an  orientation  of  the  source  such  that  its  main 
lobe  is  oriented  in  a  direction  near  the  angle  vc  should  intensify 
the  radiated  field.  For  low-loss  media  maximum  radiation  occurs 
in  the  plane  normal  to  the  axis  of  the  inclined  dipole,  such  that 
®m  -  Ve.  For  the  case  of  media  with  losses  am  will  depend  on 
these  losses. 

From  (27)  one  observes  that  maximum  radiation  occurs  for 
0=0  and  0  <  a  <  ir/2  or  for  0  =  jrandO>q>(-rr/2).For  0  =  0 
the  depencence  of  (27)  in  a  may  be  expressed  by 


A<*)  = 


sin  q  cos  q 
n2  —  1  +V^T 


(29) 


where  n  =  l/«j. 

For  q  real  the  maximum  of f{a)  occurs  for 


The  solution  of  the  integral  in  (19)  corresponding  to  the 
lateral  wave  at  the  interface  of  media  2  and  3  will  be  strongly 
attenuated  due  to  the  high  conductivity  of  the  medium  3  for 
the  cases  that  were  considered  here.  The  lateral  wave  at  the  in¬ 
terface  1  and  2  will  be  calculated  in  a  way  similar  to  that  avail¬ 
able  in  the  literature  [1]  and  [3].  Thus  the  z  component  of  the 
electric  field  will  be  given  by 

E.b  =  -Eiz  [1  +  )*“*'] 

(*  ~ni) 


tg2a„  = 


2  Re  {y^rp 
|(na  —  7)1—  1 


(30) 


This  result  is  identical  to  that  obtained  by  Staiman  and  Tamir 
[14]  for  a  half-space  geometry  since  the  contribution  of  the 
lateral  wave  in  the  forest-ground  interface  was  neglected  for  the 
practical  cases  considered  here.  This  problem  has  also  been 
studied  for  a  medium  with  three  layers  [IS]  with  the  source 
located  inside  the  intermediate  layer  and  the  observation  point 
in  the  upper  layer. 


*  S  l*>i)e2bH]m 

m=0 

where 

_  WoPz  /(*!'-  $)+*« 

Etx  = - r  e 

InkS 

b=jk2y/\-n\ 

S  =  2H  —  z  —  z'. 

Following  similar  procedure,  the  z  component  of  the  electric 


V.  RESULTS 

The  radio  loss  LB  from  an  inclined  radiating  dipole  to  a  verti¬ 
cal  receiving  antenna,  neglecting  the  influence  of  the  ground  on 
the  radiation  resistance  of  the  dipole,  is  given  by 

=  36.57  +  20  log/  +  20  log  E0  -  20  log  £*  (31) 

where  /  is  the  frequency  in  MHz,  E0  is  the  unattenuated  field 
strength  in  jiV/m  expected  from  the  transmitting  system  at  one 
mOe,  and  E^  is  the  calculated  value  of  field  strength  produced  as  a 
result  of  radiation  from  the  transmitting  sys^m  used  to  deter¬ 
mine  E0  above. 


(23) 

(24) 

(25) 

(26) 
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TABLE  I 

Characteristic  parameters  of  typical  forests 


HEIGHT  OF 

THE  TREES 

H(a) 

10 

20 

30 

PARAMETERS 

OF  THE 

e2 

1.1 

1.3 

1.3 

VEGETATION 

OjCmS/m) 

0.1 

0.3 

1.0 

PARAMETERS 

OF  THE 

£J 

20 

SO 

iO 

GROUND 

°3  («S/B) 

10 

100 

100 

Fig.  2.  Predicted  path  loss  for  a  dipole  buried  in  forest,  (a)  Tropical  “average"  forest,  t2  =  1.1,  o2  =  0.1  mS/m,f3  =  20,  o3  = 
10  mS/m,  am  =  63.4°.  (b)  Tropical  "dense"  forest,  t2  =  1.3.o2=0.3  mS/m.*3  =  50,  o3  =  100  mS/m,  om  =  46°.  (c)  Equa¬ 
torial  "dense"  forest,  t2  =  1.3,  o2  =  1  mS/m,  e3  =  50,  o3  =  100  mS/m,  am  -  25.9°.  z  =  z'  =  10  m.  Vertical  (V),  horizontal 
(H),  and  inclined  (1)  polarizations  at  6  MHz. 


Due  to  the  large  range  of  parameters  measured  for  forests  [  16]  - 
[18]  it  is  practically  impossible  to  display  curves  for  all  possible 
combination  of  these  parameters.  For  purposes  of  illustration 
three  cases  were  considered  and  are  shown  in  Table  I. 

In  Table  I  the  first  and  second  columns  correspond  to  medium 
and  dense  tropical  forests,  respectively,  while  the  third  column 
corresponds  to  values  recommended  for  the  Amazon  region 
(equatorial  dense  forest)  based  on  the  measurement  of  field 
strength  made  by  Assis  [18]  in  Tabatinga,  Amazon. 

In  Fig.  2  the  calculated  path  losses  are  shown  as  functions 
of  the  distance  between  transmitter  and  receiver  for  vertical,  hori¬ 
zontal,  and  inclined  polarizations  and  for  the  three  typical  forests 
considered  here  for  transmitting  and  receiving  antennas  located 
at  z  =  z'  =  10  m.  As  may  be  seen  in  these  figures,  the  bnsic  trans¬ 
mission  loss  is  improved  from  0.S  to  3.0  dB  by  using  an  optimum 
inclination  of  the  transmitting  dipole. 

In  Fig.  3  the  lateral  wave  powei-iadiation  patterns  for  the 
three  typical  forests  are  shown  at  a  frequency  of  6  MHz  and  for 
vertical,  horizontal,  and  inclined  polarizations.  The  heights  of  the 
transmitting  and  receiving  antennas  are  the  same  as  those  used 


in  Fig.  2.  The  enhancement  of  the  lateral  wave  and  the  increase  in 
directivity  may  be  observed  in  these  figures  as  the  polarization 
reaches  its  optimum  inclination.  This  suggests  that  the  enhance¬ 
ment  of  the  lateral  wave  may  increase  considerably  if  an  antenna 
with  a  high  gain  is  directed  along  am . 

VI.  CONCLUSION 

A  method  using  dyadic  Green’s  functions  far  the  solution  of 
the  electromagnetic  wave  propagation  in  layered  media  has  been 
developed.  This  method  presents  some  advantages  with  respect 
to  the  conventional  method  (the  Hertz  potential  method),  such 
as  1)  the  dyadic  Green’s  functions  outside  the  source  region  al¬ 
low  a  higher  degree  of  flexibility  with  respect  to  the  space  coor¬ 
dinates  and  assure  a  good  exponential  convergence;  2)  the  source 
may  contain  an  arbitrary  current  distribution  and  the  medium 
may  be  isotropic  or  anisotropic.  The  electric  fields  were  obtained 
considering  that  the  media  were  isotropic  and  homogeneous  and 
that  the  optimization  of  the  lateral  wave  depends  only  on  the 
properties  of  the  region  near  the  source.  The  propagation  of  elec¬ 
tromagnetic  waves  in  anisotropic-layered  media  has  not  been  as 
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Fig.  3.  Lateral  wave  power-radiation  patterns  Tor  dipole  buried  in  forest, 
(a),  (b),  and  (c)  are  described  in  Fig.  2.  z  =  z’  =  1 0  m.  Vertical  (o  =  90°), 
horizontal  (a  ■  0°),  and  inclined  (am)  polarizations  at  6  MHz. 
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yet  widely  explored  and  the  use  of  the  dyadic  Green’s  function 
method  may  unveil  some  interesting  results. 
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